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Proposed Silicon-based quantum computer architectures have attracted attention because of 
their promise for scalabihty and their potential for synergetically utilizing the available resources 
Eissociated with the existing infrastructure of the powerful Si technology. Quantitative understanding 
of and precise physical control over donor (e.g. Phosphorus) exchange are crucial elements in the 
physics underlying the proposed Si-based quantum computer hardware. An important potential 
problem in this context is that inter-valley interference originating from the degeneracy in the Si 
conduction band edge causes fast oscillations in donor exchange coupling, which imposes significant 
constraints on the Si quantum computer architecture. In this paper we consider the effect of external 
strain on Si donor exchange in the context of quantum computer hardware. We study donor electron 
exchange in uniaxially strained Si, since strain partially lifts the valley degeneracy in the bulk. 
In particular, we focus on the effects of donor displacements among lattice sites on the exchange 
coupling, investigating whether inter-valley interference poses less of a problem to exchange coupling 
of donors in strained Si. We show, using the Kohn-Luttinger envelope function approach, that fast 
oscillations in exchange coupling indeed disappear for donor pairs that satisfy certain conditions for 
their relative positions, while in other situations the donor exchange coupling remains oscillatory, 
with periods close to interatomic spacing. We also comment on the possible role of controlled 
external strain in the design and fabrication of Si quantum computer architecture. 

PACS numbers; 71.55.Cn, 03.67.Lx, 85.30.Vw 



I. INTRODUCTION 

Advances in quantum computing software research,0 such as the invention of factorization algorithm and quantum 
error correction codes, have pronapted active and extensive search for an appropriate physical system for implemen- 
tation of quantum computation.BI For example, there has been considerable recent interest in the study of donor 
impurities in silicon, particularly the Si:^^P system, because of the potential of the monovalent ^^P impurities to aoL 
as fundamental units of a solid state quantum computer (QC). The first Si-based QC architecture, proposed by Kane,a 
has drawn immense interest from the experimental community. In Kane's proposal the nuclear spins of phosphorus 
(■^^P) donors are the quantum bits (qubits), while the donor electrons act as shuttles between neighboring nuclear 
spins. Isotopically purifiedp^n^i host provides a quiet environment with very long coherence times for both donor 
electronic and nuclear spinaj'El since it has been known for a long timeO that the electron spin coherence time in Si is 
extremely long. Two-qubit opecations, which are required for a universal quantum computer, involve precise control 
over electron-electron exchangeQ and electron-nucleus hyperfine interactions. A closely related alternative design was 
proposed later ,□ involving ^^P donors in Si/Ge heterostructures, with the spins of the donor electrons (rather than the 
nuclear spins) serving as qubits. The basic ingredients underlying the Si QC proposals are schematicaUy presented and 
described in Fig. ^ There are also other recent QC proposals based on Si or Si/Ge host materials .ffl Although these 
are interesting and promising proposals, particularly due to the extensive infrastructural base available in the existing 
Si microelectronics technology, there are formidable experimental challenges to be overcome in device fabrication, 
coherent control, system integration, and single spin measurement, among other important issues, before a Si-based 
QC architecture could become an experimental reality even on a laboratory scale. 

Recently, there has been significant progress in the fabrication of donor arrays in silicon single crystal, using both-a. 
"top-down" approach with ion-implan±ation, and a "bottom-up" approach with MBE growth and STM technique.EHl 
Theoretically, we have recently shownliil that electron exchange in bulk silicon has fast atomic scale spatial oscillations 
due to the valley interference arising from the very special six-fold degeneracy of the bulk Si conduction band. These 
oscillations place heavy burdens on the device fabrication and coherent control because of the extremely high precision 
requirement for pJiacing each donor inside the Si unit cell and/or for controlling the external gate voltages. 

Several authoraJ-El have pointed out that Si-based quantum computing architectures involving epitaxial_Si/Sii_2;Gea; 
heterostructures possess some practical advantages over the original scheme based on relaxed bulk SLd One of the 
motivations is the need for a barrier material separating the Si host containing "^^P donors from the conducting gates 
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(see Fig The Si/barrier-material interface must be free of structural defects. Although high quality Si/oxide 
barriers can be fabricated, Si/Sii_a;Gea; and Si/Sii-^Ca; are strong candidates for QC architectures because of the 
good interface quality associated with their epitaxial nature. These alloys have advantages over the oxide barriers (e.g. 
8102) in terms of chemical compatibility among the group IV elements, which share the diamond lattice structure, 
while still providing sufficiently high barriers for the donor electrons to prevent leak into/from the gates. In spite of 
sharing the same lattice structure, C, Si and Ge present very large lattice mismatches (lattice parameters 3.57 A, 5.43 
A, and 5.66 A, respectively), leading to the presence of strained layers in heterostructures containing commensurately 
grown materials. The highly strained nature of these epitaxial systems (arising from their large lattice mismatch) in 
fact is crucial in producing high quality defect free interface since dislocations are avoided due to the strain build-up. 

Given our earlier finding of band-degeneracy-induced strong oscillations in the donor exchange energy (as a function 
of individual donor positioning within the Si unit cell) of bulk Si, it is natural to ask how donor exchange behaves in 
strained Si heterostructure or quantum well systems, where one must now incorporate strain effects in the exchange 
calculation. In this paper we carry out such a calculation of strain effects on exchange energy, and emphasize an 
additional feature of '^^P-doped strained Si quantum wells (QWs) favoring such architectures over architectures based 
on relaxed bulk Si for QC implementation. Recall that two-qubit operations of a Si QC rely on the exchange coupling 
between the wave functions of neighboring donor atoms. The six-fold degeneracy in the conduction band edge of Si 
leads to oscillations in the exchange couplingjlil which pose formidable experimental problems in controlling both the 
fabrication (in terms of donor positioning within the crystal) and the actual operation (in terms of surface gate bias) 
of the QC. We demonstrate below that these problems are partially solved for impurity pairs in strained Si grown 
over Si/Ge alloys due to altered lattice and band structure. Unfortunately, however, we find that some parts of the 
exchange oscillation problem persist even in these strained systems as we discuss in this paper. 

In the following, we first discuss how to describe theoretically the effects of strain in silicon. We then discuss the 
calculation of the donor exchange energy via an effective Heitler-London approximation. Finally we present our results 
and discuss their implications and validity. The details of our Heitler-London calculation are given in the Appendix. 

II. THEORETICAL APPROACH 
A. Strain in Si layers 

Epitaxial growth techniques such as molecular beam epitaxy (MBE) allow the fabrication of lattice-mismatched 
heterostructures |fcee of misfit defect generation when the layers are sufficiently thin (and the lattice mismatch is not 
extremely large) .113 The mismatch is completely accommodated by uniform lattice strain and the interatomic distances 
parallel to the interfacial plane, i.e. the "effective" lattice constants in the plane, remain equal to the equilibrium 
value of the substrate material. The lattice constants perpendicular to the interface adjust appropriately in order 
to minimize the total elastic energy in the strained layer. This so-called commensurate growth lowers the energy 
of the interfacial atoms at-the expense of increasing the elastic (strain) energy stored in the chemical bonds in the 
coherently strained layers.li3 We refer to these layers as uniaxially strained. Depending on whether the interatomic 
distances parallel to the interface are larger or smaller than the relaxed equilibrium value, the strain is called tensile 
or compressive. The structural variations involved in this process modify the electronic properties of the layered 
materials and consequently the substitutional impurities in the strained layers. p 

Here let us consider the important example of Si grown over a relaxed Sii-^^Ge^: suJasjtxate on a (001) interface.cl 
Successful commensurate growth of such heterostructures has been previously reported.Ejil3 The average alloy lattice 
parameter, which follows Vegard's law very closely,t3 defines the Si layer lattice parameter parallel to the interface: 
aj| = (1 — a;) asi -I- x acc- The lattice parameter perpendicular to the (001) interface is also modified in the strained 
Si layer with respect to its equilibrium value. According to macroscopic elasticity theory, aj_ = asi -I- 2(ci2/cii) (agi 
- suA. Here C12 = 16.577 dyn/cm^ and cn ~ 6.393 dyn/cm^ are elastic constants with values given for relaxed bulk 
Si.llll The elasticity theory predictions for the crystal lattice structures in Si/Ge heterostructures were confirmed by 
ab initio calculations by Van de Walle and MartinEj For an alloy with x = 0.2, the above relations yield ay = 5.474 
A and a_L = 5.396 A, i.e., distortions of less than 1% (-1-0.8% and —0.6%) with respect to the equilibrium value of 
asi = 5.43 A. 

Uniaxially strained Si has lower symmetry than relaxed Si, and thus a different reciprocal lattice (body centered 
tetragonal instead of body centered cubic) and Brillouin zone (BZ). Herring and VogttSl have shown that the small 
tetragonal distortions cause changes in the conduction band minima of Si that can be quantitatively described by 
shifts in the energies of the local minima of the six valleys, while the reciprocal-space positions and shapes of the 
constant-energy surfaces remain unchanged to first order in strain. The energy shifts are given in terms of the amount 
of distortion and of the uniaxial and dilation deformation potentials, and S^, respectively. For Si grown over 



2 



relaxed Sii-^Ge^; (001) substrate, leads to uniform shifts of the six valley minima. Therefore, as discussed below, 
only is relevant for our study of '^^P donor wave functions. 



B. Effective Hamiltonian for a donor electron in strained Si 



The conduction band edge of relaxed bulk Si consists of six degenerate minima located along the (100) directions, 
at about 85% awayriom the BZ center, towards the zone boundary at X points. Within the Kohn-Luttinger envelope 
function approachjij the ground state for a donor electron in Si is also six-fold degenerate, and the corresponding 
states are labeled by the reciprocal-space points at which the valley minima occur. For definitenessT,we follow the 
sequence: /x = 1, 2, . . . , 6 <-^- z, — z, y, — y, cc, — x, with the envelope functions given by (e.g. for fj, = ±z)i3 

F^^ir) = -L= ,-[(-^+.^)/a^+^V^.Y/^ . (1) 

The effective Bohr radii for Si from a variational calculation are a — 25.09 A and 6 = 14.43 A.0 

A substitutional impurity breaks the translational symmetry of the host crystal, leading to intervalley scattering 
effects known as the valley-orbit interactions,E3'Ell which split the unperturbed six-fold degenerate donor electron 
ground state into a singlet (^i), a triplet (Ti) and a doublet (_£), with th*, ground state being a non-degenerate 
state of Ai symmetry. We follow the perturbative treatment of Baldereschi,E3 which starts from the six degenerate 
ground-state functions calculated within the one-valley Kohn-Luttinger variational approach. The low energy donor 
electron wave functions can be expanded in the basis spanned by the six valleys 
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V'(r) = ^a^F^(r)</)^(r), (2) 

where (j)fi{v) = M^(r)e*'''' '" are the pertinent Bloch wave functions and = 1- Each coefficient is in 

general referred to as ^l- valley population^ while a set of 6 coefficients defines the so-called valley composition of a 
donor electron state. 

To determine the donor electron wave function, we include as perturbations two types of intervalley coupling, for 
valleys on perpendicular symmetry directions (e.g. x, z) and for those on opposite symmetry clirections (e.g. z, 
—z). We represent these couplings by — Ac and — Ac(l -I- 6), respectively. The experimental valuea£3 for the relative 
splittings among the Ai, Ti and E levels for ■^^P donors in Si are obtained as Ac = 2.16 meV and S = —0.3. 
Taking valley-orbit scattering into consideration through these parameters leads to a binding energy of 40 meV for 
the Ai ground state, in quite a good agreement with the experimental value of 45.5 meVJSS More accurate estimates 
of the binding energy were obtained through non-perturbative variational treatments,Ej which however rely on a 
spherical-band approximation and adopt s-like hydrogcnic trial envelope functions. Since our goal here is to obtain 
a realistic description of the ground state wave function for the donor electron rather than the precise values of the 
binding energy, we employ the perturbative approach with anisotropic envelope functions. For our considerations of Si 
QC hardware architecture, which focus on electron exchange due to wavefunction overlap, the perturbative envelope 
function approximation that we employ is quite adequate. 

Uniaxial strain along [001] induces different valley shifts for the two local minima along the z axis (/i — 1,2) 
compared to the other four along x and y axis (/i — 3,4,5,6). Within the subspace spanned by the electron wave 
functions at the six conduction band minima {.^^(r) (/)^(r)}^=i,6, the donor ground and lower excited states may be 
conveniently obtained as the eigenvectors of the effective perturbation Hamiltoniaa in the form of a 6 x 6 matrix 
which, for ''^P donors in a Si QW uniaxially strained along the z axis, is written afj 

/ 1+6 1 1 1 1 \ 

1+S 1111 
1 1 Q+S 1 1 
1 1 1+5 1 1 
11110 1+5 
V 1 1 1 1 1+5 J 



H — Hi,n "~h -^strain 
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Note that the order of the rows and columns of the perturbation Hamiltonian H foUows the sequence in /i given above 
Eq. (|l|). Furthermore, since our concern here is the donor wave functions, any shift in H proportional to the identity 
has been neglected, so that the eigenvalues of H give the correct relative splittings among the eigenstates, but not 
the binding energies. The first term of the perturbation Hamiltonian, H^o, gives the valley-orbit scattering due to the 
presence of ^^P donors in unstrained Si discussed aboverJ The second term, iJstrain, introduces the relative energy 
shifts due to uniaxial strain along [001] in terms of a dimensionlessi-scalar valley strain parameter x defined in terms 
of strain parameters and the valley-orbit scattering matrix elements 
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Here the uniaxial strain parameter :i„L3 is approximately 8.6 eV for the Si conduction band edge, while ct — — 
Cxx is the difference between the diagonal components of the strain tensor for uniaxial strain along z: — a_L /asi-1 
and Bxx = ay/asi — 1, so that the valley strain parameter x is linearly related to the Sii-^Ge^: alloy composition x as 
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Since S„ is three orders of magnitude larger than Ac (with S„ ~ 10 eV and Ac ~ 2 meV), relatively small changes in 
X may lead to important shifts in the energies of the valleys. For example, for Si grown over a relaxed Sii_a;Gea; alloy 
with X — 0.2 on a (001) interface, lattice distortions (smaller than 1%) lead to a valley strain parameter x ^ —20. 
For Sii-ajGca; alloys of arbitrary composition, we estimate xi^) — —95a;. Negative values of x indicate tensile strain 
(a|| > asi), favoring the z valleys energetically. Commensurate growth of Si over Sii-^^Ca, or SiC, which is more 
difhcult experimentally due to the extreme lattice mismatch between the materials, would lead to compressive strain 
(ay < asi), so that x > 0. The lattice mismatch between C (diamond) and Si is 34% (3.57 A versus 5.43 A), as 
compared to 4.2% for Ge with respect to Si (5.66 A versus 5.43 A). When we replace acc by ac in Eq.(^), the 
estimated valley strain parameter for Si grown over relaxed Sii-^^Cj, would be x ~ 800a;. This is actually ja-, lower 
bound for x, since this alloy's average lattice parameter is somewhat smaller than predicted by Vegard's lawJl3 

A symmetry argument based on the differentiated fit of the six envelope functions to the Si host geometry and 
boundary conditions may be used to include the effect of the confinement potentials due to the alloy regions in 
Sii_a;Ge3;/Si/Sii_2,Ge3; heterostructures, or due to an interface with any potential barrier (e.g. Si/Si02) within a 
few effective Bohr radii from the impurity. Assuming that all Si/barrier interfaces are perpendicular to the growth 
direction z, the components with F±z envelopes are favored energetically, regardless of the value or sign of x, since the 
smaller effective Bohr radius b along z guarantees a less significant penetration of the wave function into the barrier 
regions as compared to the other envelopes with the larger Bohr radius a along z. This effect is phenomenologically 
included in the third term of the donor electron Hamiltonian H through an energy shift A. The parameter A 
is alwaySpBOsitive, and its value depends on the barrier height and on the impurity position with respect to the 
interface.^ We estimate it here to be up to ~ 10 meV. 

In addition to valley energy shifts, the presence "ptan interface perpendicular to z also leads to coupling between the 
±z valleys. Surface-induced intervalley scattering is introduced phenomenologically in through the pacajsieter 
il. This coupling is in general complex, and leads to energy shifts smaller than 1 meV in realistic situations.EZrEa 



C. Single donor electron states and energy spectrum 



In the absence of strain and other perturbations, the single donor electron Hamiltonian consists only of the valley- 
orbit coupling, i.e. H — Hyo- If (5 = 0, so that all the valley-orbit scattering matrix elements are the same, the 
eigenstates of H^o are a singlet with energy — 5Ac and a quintuplet with energy Ac. When S 0, the quintuplet 
splits into a triplet and a doublet, so that a single donor electron has the state sequence of a singlet, a triplet, and 
a doublet, which have energies —(5 -I- S)Ac, (1 -I- S)Ac, and (1 — (5)Ac, respectively. For Si, where 6 ~ —0.3 and 
Ac ~ 2.16 meV, the singlet is at — 4.7Ac, the triplet at 0.7Ac, and the doublet at 1.3Ac. The first excited state here 
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is a triplet that is 5.4Ac ~ 11.7 meV above the ground singlet state. The valley compositions of these states show that 
the ground singlet has Ai symmetry, corresponding to a symmetric superposition of all six valleys: ■^(Ij Ij li Ij Ij l)j 
while the triplet and doublet have Ti and E symmetries as expected.ll2l 

When the Si lattice is uniaxially strained along the z axis, H = Hy^ + Hstmin- Now the degeneracy in the triplet and 
doublet states are also lifted. For x < 0, the ground state energy is Ac —{2 + 6) + x/2 — (3/2)-\/x^ + 4x/3 + 4 : It 
is an admixture of the Ai state with one of the components of Ti, with valley composition (ai, ai, 013, as, as, as) and 
as = ai/ y'x^ + 4x/3 + 4 — (x + 2/3) . Notice that the ground state approaches (1, 1, 0, 0, 0, 0)/-\/2 as x —oo- In 
this limit , the first excited state is one of the original triplet states (in the absence of strain) , with energy Ac(l + <^+2x) 
and valley composition (1, —1, 0, 0, 0, 0)/-\/2. The energy splitting between these lowest states approaches 2{l+S)Ac ~ 
lAAc ~ 3.02 meV. For x > 0, the lowest energy states are admixture of Ai with the other two components of Ti: 
The (nondegenerate) ground state has energy Ac h (2 + <5) ~ x/2 - (3/2) Vx2 + 4x/3 + 4] ^ Ac (-3 - ,5 - x) , while the 
first excited state (actually a doublet) has energy Ac(l + ^ — x)j with a splitting from the ground state approaching 
(4 + 26)Ac ~ 3.4Ac ~ 7.34 meV as x ^ +00. Hydrostatic (dilation) strain causes a rigid energy shift in the 
six valleys, keeping the splittings and valley compositions of all eigenstates unchanged. Since our focus here is to 
understand wavefunction overlaps, we do not consider effects of hydrostatic strain. 

In Fig. panel (a), we plot the energy splitting between the first excited state and the ground state of a single donor 
electron as a function of the valley strain x- This quantity is relevant in estimating time scales which determine the 
adiabatic condition in time-dependent processes driven by the gate potential variation. As we have shown above, the 
energy splitting becomes smaller when the bulk ground and excited states mix to form strained ground and excited 
states. The variation of these energies can be seen in Fig. 1 of Ref. |2^. It is also clear from that figure why the curve 
we plot here has a cusp at x = 0: A level crossing occurs at x = for the first excited state. 

The strength of the hyperfine coupling between the donor electran and nuclear spins is also important for the Si QC, 
since it is invoked in both single qubit and two-qubit operations.u In terms of the donor electron state composition, 
it can be expressed as 
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A{x)/Ao = I 



E 



< 1 , (6) 



where Aq is the value of A{x — 0). This is a differentiable function of Xi so its peak at x = is smooth, as is shown 
in panel (b) of Fig. ||. Furthermore, for x —00, the ground state composition approaches l/\/2(l, 1, 0, 0, 0, 0), so 
that A{—oo)/Aq = 1/3; while for x 00, the ground state approaches 1/2(0,0, 1,1, 1, 1), and A{oo)/Ao = 2/3. The 
curve in Fig. ||(b) clearly approaches these two limits. Results summarized in Fig. || indicate that donors in uniaxially 
strained Si will have "inferior" single-electron properties compared to the relaxed material for the purpose of QC 
operations, requiring longer times to perform an adiabatic evolution of the electron state and providing a reduced 
hyperfine coupling. 

When the z-direction confinement is added, it introduces very small energy shifts and additional off-diagonal 
coupling between z and — z directions. In the presence of Hyo and/or iJstiain, the effects of Hz are negligible for all 
relevant properties discussed here. Thus we do not include this term in our calculation below. 

We emphasize the following general properties of the spectrum of H in Eq. (^) when Ac > 0: (i) the ground 
state is always non-degenerate, and (ii) the ground state valley populations in opposite symmetry directions, e.g. 
fi = 1, ly = 2, always satisfy la^^P = Icki^p. 



D. Heitler-London approximation 

In this paper we use Heitler-London (HL) approximation to calCiWilate the inter-donor exchange coupling as the 
energy difference between the two-electron singlet and triplet statesEj Considering a pair of donors at TIa = and 
Rs = R and assuming i? = |R| 3> a,6 (the effective Bohr radii), the energies of the singlet and triplet states are 

where Eq is the single impurity electronic energy, while Hq, Hi, and S are integrals involving two electronic wave 
functions of the form (H) centered at the origin and R, and are given in the Appendix. Contrary to the case of the H2 
molecule, where these integrals are isotropic and smoothly decaying functions of the relative distance R, in the case 
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of Si they are extremely anisotropic and sensitive functions of R.0 The HL expression for the exchange splitting in 
Si is 

J(R) =Et-Es 

= E [ E l^nc^c'Pe^^K-'''^-''] lapPla.P Jm.(R) cos(k^ - k,) • R, (8) 

K,K' 

which is derived in the Appendix, where the explicit expression for J^^v{^) is also given. The second summation 
(the part within the square bracket) in (||) refers to the reciprocal lattice expansion of the periodic part of the Bloch 
function, w^i(r) = c^e*^ "", and is identically unity when R is an fee crystal lattice vector. 

A prominent feature of the above expression for the exchange coupling J(R) is the presence of the fast oscillatory 
terms cos(k^ — kj^) • R, which are periodic with wavelengths of the order of the atomic spacing (^ 5A in Si). The 
coefficients J^i, (R) are slowly varying functions of R as they are integrals containing the envelope functions (R) , 
whose characteristic decay length is the effective Bohr radius of the order of 20 A. 



III. RESULTS 



We have performed HL calculations of the donor electron exchange for different configurations of a pair of donors. 
Results (Figs. ^^) presented here are obtained by taking the summation over reciprocal lattice vectors to be unity: 
K' kKpl''K' P^*^^^^ ■* ^ = 1' i-^- term within the square bracket in Eq. (||) is taken to be unity which is 
equivalent to calculating J(R) only for those values of R which correspond precisely to R being an fee lattice vector. 
The continuous lines shown in the figures correspond to the free-electron ilifnit « ^k,o, which is a first order 
approximation to the lowest conduction band in homopolar semiconductors. e3 We investigate the exchange coupling 
in uniaxially strained Si as we vary the inter-donor distance, the inter-donor direction relative to the strain axis 
(defined to be the z axis), and the sign and the magnitude of strain itself. Results presented in Figs. ^ to |^ give a 
comprehensive account of the general trends and effects to be expected, providing new qualitative and quantitative 
insight into the problem of donor electron exchange coupling in multivalley semiconductors. 

In Fig. II we plot the exchange coupling J as a function of the inter-donor distance R and valley strain x for a 
pair of donors along the [100] axis. Notice that as x increases from negative to positive values, the exchange coupling 
becomes increasingly oscillatory. Such behavior can be easily understood from the expression of J(R) in Eq. (H) and 
from the donor ground state valley composition. For x ^ with strain along the z axis, the dominant components 
of the ground state are 

V'(r) ^ aiF.e^'^-'- + a^F^^e-^'^'-' , (9) 

where ax = a2 ^ l/\/2, so that the largest oscillatory contribution to J(R) comes from terms with factors cos(fcz ± 
k^z)Rz with k^z = —kz- The only effect of varying the x-component of R is to vary the coefficients J^^u, which are 
slowly decaying functions of |R|. As x increases toward zero, the contribution from the four components as, . . . , ag 
grows, and fast oscillations are superimposed on the overall decaying envelope of J{Rx)- For x ^ 0, these four 
components dominate the ground state valley composition, and the oscillations arc enhanced with respect to relaxed 
bulk Si. 

From the perspective of predicting and controlling the donor pair exchange coupling, results in Fig. || indicate that 
the optimal degree of strain is in the range x ~ since this leads to a simple hydrogenic-type behavior for the 
particular type,pf relative positions (along [100] axis) of a pair of donors, as was presumed in the original Si and Si/Ge 
QC proposals.H'Q Moreover, the value of J is considerably enhanced (by about 100%) with respect to values in the 
relaxed bulk Si, a feature with considerable practical advantages in terms of increasing the speed of gate operations 
in the QC. 

To study cases when the two donors are not positioned along the [100] axis, we first choose a reference configuration, 
with a donor pair along [100] direction. We choose to work with uniaxially strained Si with x — ~20 from now on as 
it corresponds to the realistic situation of Si grown over relaxed Sio.8Geo.2- We arbitrarily pick the relative position 
Rx ~ llOA marked by the solid triangle in Fig. ^ with respect to which various displacements are considered for 
one of the two donors. In Fig. ^ we plot the exchange coupling J as one of the donors is displaced from the reference 
site along either the y ([010]) or the z ([001]) axis. The upper curve, corresponding to displacements along the y axis, 
reflects essentially the same qualitative behavior as in Fig. ^ for x = —20. However, when the relative displacement 
R acquires any nonzero Rz component, exchange becomes rapidly oscillatory, as indicated by the triangles in the 
lower curve. The same trends are seen in Fig. S, where displacements along the [110], [Oil] and [101] directions are 
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considered. Among those, only displacements along the [110] diagonal direction lead to a monotonous variation of J 
with distance without any oscillations. The other two types of displacement involve a nonzero Rz component in the 
relative position R, leading to fast oscillations in the exchange energy. 

The overall effect of uniaxial strain in the exchange coupling of a donor pair in the x-y plane is summarized in 
Fig. ^. There we consider all possible in-plane substitutional donor pairs with intcr-donor distances in the range of 
90 A < R < 180 A. For impurity pairs in uniaxially strained Si (x = —20), the exchange is a single- valued function 
of R, decaying smoothly with distance, while for unstrained Si (x — 0) the exchange is significantly reduced with 
respect to the strained results, and is sensitively dependent on the relative position R. 

Notice that in all the results presented here the oscillatory behavior of J is almost always accompanied by partial 
cancellations, meaning that in the oscillatory regime, the values of J are always smaller than in the corresponding 
non-oscillatory situations. This is clearly evident in Fig. ^ comparing curves with increasing values of x; or in Figs. | 
and ^, comparing curves with zero and nonzero Rz components for the same inter-donor distance |R|. This result 
can be mathematically confirmed from expression (^, where the cos(k^ — k,y) • R factors are responsible both for 
oscillations and for partial cancellations [j7'/i,y(R) > for arbitrary fj, and i^]. It is somewhat different (and perhaps 
counterintuitive) from an ordinary two-wave interference. One aspect of our results presented in the figures (Figs, y- 
^ that needs to be emphasized is that our exact HL calculations are done only for discrete values of the relative 
displacement or inter-donor distance, and the smooth curves are guides to the eye. As mentioned in the beginning of 
this section, the smooth continuous curves in our results correspond to an effective free electron approximation. 



IV. SUMMARY AND DISCUSSIONS 



A. Approximations 

We now briefly review the main approximations involved in the theoretical approach adopted in this paper, and 
discuss various possible improvements. In our perturbation theory calculation, the only effect of strain is to change 
the valley populations of the donor electron states through the diagonal elements of the six-valley Hamiltonian H. 
We do not consider strain-dependent changes in the effective masses, the inter-valley scattering matrix elements Ac 
and (5, or the BZ locations for the conduction band minima. Since the lattice parameter changes are quite small — less 
than 1% as we mentioned in the beginning of Sect. II, we believe that our approach is well-suited for the modestly 
strained systems considered here. For highly strained systems (such as Si grown on Sii_a;Ca; with high concentration 
of carbon), the separation of H into Hyo and i?stiain may not be justifiable, and all matrix elements should be 
consistently assessed to determine how they depend on strain. 

In our envelope function approach, we have adopted the variational form and parameters obtained from the one- 
valley Kohn-Luttinger treatment. Calculations incorporating the valley-orbit interactions and strain in the original 
Hamiltonian, and treating the effective Bohr radii a and b and the valley compositions as independent variational 
parameters (imposing symmetry and normalization constrains), would have led to a better quantitative description 
of the ground state binding energy and possibly the wave function. Such calculations for donors in .Si have been 
performed to determine the binding energies using only spherically symmetric trial wave functions,LirEj which is not 
adequate for the present study. Our interest in this paper is not in obtaining extremely accurate donor binding 
energies, but in calculating reasonably good wavefunctions so as to obtain reliable donor exchange energies. 

Atomistic ab initio studies based on density-functional schemes can overcome the limita±ipns of the envelopiei 
function/effective- mass theory and have been employed in the literature to study effects of strair£3 and alloy disorder .113 
The difficulty in describing a single shallow impurity originates from the underlying supercell approach with periodic 
boundary conditions, which leads to spurious effects due to interactions between the impurity and its periodic replicag. 
Extrapolation schemes towards an infinite supercell have been proposed for the single impurity binding energy,E2l 
although periodic boundary conditions cause additional difficulties for a reliable description of the asymptotic behavior 
of the wave functions. 

The HL approximation adopted inj-this paper is the simplest, approach to calculate exchauge coupling in a variety 
of situations, such as atomic lattices,t3 random donor arraysj^j and double quantum dots.L3'Lj Improvement can be 
made through various ways. For instance, it is quite common to use the valence bond theory (also called the molecular 
orbital theory) or the configuration interaction (often combined with a self-consistent mean field calculation) approach 
to obtain pinre accurate spectra (and thus exchange interaction when specifically pursued) for molecules or quantum 
dots.E3Hi2a For donor arrays in Si (relaxed or strained), such calculations will be complicated by the Si band 
structure. However, as we have shown in the current study, as long as the effective Bohr radius is much larger than 
the lattice spacing, one can separate the fast oscillating (due to the Bloch functions) parts of integrals for the matrix 
elements from the slowly varying envelope function parts, thus simplifying the calculation. Indeed, there should be 
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no significant theoretical difficulty for a more complete molecular orbital calculation for donors in Si. In the current 
study, our focus is on the qualitative behavior of the exchange coupling between two donors, therefore we have not 
attempted to improve upon the HL approximation. This can be easily done following earlier calculations in GaAs 
quantum dotsu if the experimental development warrants such improvements of the theory. In addition, the proposed 
Si QC jsipuld operate under a uniform magnetic field of 1 to 2 Tesla applied along the z axis (see Fig. Previous 
studiesH indicate that uniform magnetic fields of up to 2 Tesla have no qualitative effect in the exchange coupling 
of electrons in double quantum dot bound states. For a pair of donors, a reduction in J due to the shrinking of the 
envelopes perpendicular to the magnetic field direction and magnetic phase interference should be expected, but this 
should not lead to any qualitatively different effect. 

As we have mentioned before, in our calculation we have approximated the Bloch functions by the-corresponding 
plane wave functions, which is in essence a free electron approximation for the conduction band of Si.EJ The approx- 
imation here amounts to keeping only the cq term (and let cq = 1) in the reciprocal lattice expansion of the periodic 
part of the Bloch functions. The plane wave approximation actually obtains the exact results for lattice points on 
the same fee lattice, where e*'^^^^ = 1 so that the sums over reciprocal lattice vectors K are normalized to 1. 
However, for an arbitrary R, our calculation should mostly give an upper bound to an exact HL calculation. If the ck 
coefficients have a broad distribution, the sums would tend to be much smaller since they involve the fourth power of 
Ck as is shown in Eq. (H). Thus the continuous lines in our figures (Fig. ^ to Fig. 0) should be taken more as guides 
to the eye than as the actual values of exchange, which can be much smaller. As mentioned before, we carry out 
our exact Heitler-London exchange coupling calculations only for discrete values of inter-donor separation (triangles, 
squares, etc. in the figures) whereas the continuous curves represent the free electron HL approximation. Further 
calculations are underway to go beyond the free electron approximation to address this issue, and those results will 
be presented elsewhere. 



B. Valley degeneracy effects 

Our results for Si uniaxially strained along the z axis indicate that, for x ^ ^5, inaccuracies in the positioning 
of substitutional ^^P donor pairs within the x-y plane produce only small effects in the electron exchange coupling. 
Naturally, the next question is whether it is possible to avoid the oscillatory behavior due to lattice-parameter-scale 
inaccuracies in the component of the relative position between a pair of donors. As noted in Sec. HI, these 
oscillations in J(R) are due to the specific form of the ground state wave function (|^), which dictates that the 
contribution from the two valley minima are equally important. This result is in turn a consequence of the degeneracy 
of the and — valleys in bulk Si, whose energies are lowered equally with respect to the four perpendicular 
valleys by uniaxial tensile strain. The same effect occurs-iin the two-dimensional electron gas (2DEG) in (001) 
Si/Si02 inversion layersEj and Si/Sii-j^Gcj: quantum wellsEj where, due to the confining potential at the interfsj:]©- 
the lowest electronic^ubbands originate from the bulk Si k^ and — k^ valleys alone. It has been shown theoreticallyE3'c2l 
and experimentallyE3c3 that this valley degeneracy may be lifted in heterostructure-type configurations, leading to 
experimental valley splittings of the order of 0.1 meV. Such splittings refer to the unperturbed doubly-degenerate 
ground states, and the reailting singlet ground state is not one of the k^ and — kx valleys, but an equally weighted 
admixture of these states£3 In the present case, Hyo alone leads to a nondegenerate ground state, which is 12 meV 
below the first excited state. Under the ifstrain perturbation, the ground state remains nondegenerate, with an energy 
separation to the first excited state reduced, but still larger than 3 meV. Further perturbations, e.g. that produces 
the valley splitting observed in the 2DEG, have negligible quantitative impact in any of the results obtained here. 
Moreover, perturbations of the form always lead to an equally weighted admixture of the k^ and — k^ valleys, 
thus donor electron valley composition always has contributions from these two valleys with the same weight, i.e. 
|ai| = \a2\ in (^. 

We note that for Ge, where four conduction band minima occur at the inequivalent L points in the BZ, (111) tensile 
strain leads to a nondegenerate band minimum. Since the lattice parameter of Ge is larger than that of Si and CL 
this possibility within group-IV semiconductors remains restricted to Ge-rich strained alloys grown over relaxed Ge,u 
with the additional complication of bringing alloy disorder effects into the non-degenerate ground state layer where 
the donors are placed. We do not consider Ge in this paper since there is no particular experimental effort directed 
toward the design of a Ge-based QC. 
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C. Estimates for static electric field effects 



Static electric fields generated by the top gates can be used to slightly modify the donor electron wavefunction in 
its shape and weight distribution, so that such fields can be used to overcome lattice-parameter-scale inaccuracies 
in the Rz component of the relative position of a pair of substitutional impurities. A crude estimate of the-uniform 
static electric field effects in the donor bound states may be inferred from results for the hydrogen atom.eS To the 
first order in the field, the H ground state wave function in the presence of a uniform electric field £ along z — r cos 9 
is given by 



V'H(r) = 



1 



-r/ao 



£ f r2 
1 - - I aor + y I cos ( 



(10) 



where aq = 0.529A is the Bohr radius, and e is the electron charge. The expectation value of the electron position 
is thus shifted by the field from (r) = to (r)H = — |aQ£/e, which corresponds to an average displacement of 
4.63 X 10"^ A for £ = 10 kV/cm. This negligible shift is greatly enhanced for donor bound states in Si. Assuming a 
spherically averaged envelope with an effective Bohr radius a* = {o^bY^^, we get (r)gi — e(a*/ao)'^(r)H — 7.4xlO^(r)H, 
where e = 12.1 it the Si dielectric constant. An average position shift of one monolayer (IML = asi/2 = 2.7 A in Si) 
thus corresponds to £ = 8 kV/cm. In the current context, it is plausible to expect that a uniform electric field of tens 
of kV/cm may compensate for atomic misplacements of a few monolayers along the field direction. 

It is interesting to note that, according to Eq.(p^), the applied field would have no first-order effect on the electron- 
nucleus hyperfine coupling, which is proportional to the electronic density at the nuclear site, A ~ |V'(0)p. Since 
the A-gates are primarily meant to reduce this coupling by pulling the electron wavefunction away from the donor 
towards the barrier, fields much stronger than those quoted above are required to perform such operations. In the. 
strong field limit, donor ionization starts to take place. The field dependence of the dissociation rate for atomic H isE^ 



AN 



^1/4 



exp 



2e 
'3£al 



(11) 



which means that dissociation becomes appreciable for field values in the order of £h ~ e/ap = 5 x 10^ kV/cm. Under 
the spherically averaged envelope assumptions above, we estimate that donor ionization in Si would require fields of 
the order of £si = [{ao/a*)^ /e]£}i = 5.3 x 10~^£h ^ 300 kV/cm, which is one order of magnitude larger than the fields 
required to compensate for atomic displacements of a few ML. 



D. Fabrication-related aspects 

We have shown in this paper that strained Si grown commcnsurately over relaxed Sii^ajGe^; (001) alloys presents 
clear theoretical advantages over relaxed Si as a host material for "^^P donor impurities in Si-based QC architectures. 
Since the same is not true for strained Si grown on Sii_a;Cj: substrates, here we focus on growth of strained Si over 
relaxed Sii_3;Gea; alloys. One important limiting factor in the fabrication of strained Si layers is the so-called critical 
thickness, he, the thickness above which the strained layer relaxes by forming topological defects, typically misfit 
dislocations which oriaujate at the interface and propagate into the bulk.Ej Theoretical estimates through an energy 
minimization criterionlij show a sharp decay of he with the alloy Ge content, from over 1000 A for x = 0.05 to 100 
A for X = 0.35. Recent advances in epitaxial growth of Si/Ge heterostructures have led to the fabrication of samples 
of 100 A thick Si layers over relaxed Sii-ajGe^; alloys with x up to 0.3jl3 with interfaces of exjteemely good quality 
and free of dislocations. Also, 150 A thick layers over x = 0.2 alloys were recently produced.ta In both cases, the 
layer thickness is below (but very close to) the theoretical upper bound for he- 

Substrate layers with relatively high Ge-content {x ~ 0.2) serve two useful purposes in the context of QC architec- 
ture: (i) the strained Si layer leads to improved behavior of exchange coupling between donor impurity states; ^ii) the 
relaxed alloy layer serves as a barrier separating the Si host containing ^^P donors from the conducting gates.Q Such 
values of alloy composition x, leading to valley strain x ^ ~20, also impose some limitations in the implementation of 
the Si-based QC. The magnitude of strain means an upper bound of around 200 A for the thickness of the strained Si 
layer, a reduction in the hyperfine coupling, and a reduction in the energy splitting between the ground and the first 
excited state by a factor of 3 with respect to the relaxed Si architectures. Alloys with lower Ge content {x ~ 0.05) 
with valley strain parameter x ^ ~5 provide comparable electronic performance (actually somewhat improved perfor- 
mance in terms of the strength of the hyperfine coupling) and greatly enhanced values of he, over 1000 A. This would 
obviously facilitate most aspects of the design and fabrication of the Si layers. On the other hand, the low-Ge-content 
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alloys do not provide a sufRciently high potential barrier, thus one would have to use the more conventional Si02 
barriers between the strained Si and the top gates. 

Our results refer to structures of strained Si commensurately grown over relaxed random alloys. The inevitable 
substrate alloy disorder can pose some potential difficulties. Indeed, some randomness in the lattice spacing and thus 
in strain is bound to occur at the alloy/Si interface. However, this effect quickly decays inside the Si layer. Since we 
consider donors located as far as 100 A (~ twenty lattice spacings in Si) away from the interface, any alloy disorder 
effect should be negligible, leading to a uniformly strained Si host for all practical purposes. On the other hand, if a 
device is based on an heterostructure predominantly formed by Sii-^Geaj/Sii-^Gej, strained layers (alloy on alloy), as 
considered in Ref. ^, and the donors are located inside an alloy layer (see Fig. ^b)), random strain and local chemical 
disorder would presumably play a much more important role in determining the donor electronic properties, including 
exchange and hyperfine couplings. The electronic properties of a single impurity is sensitive to local composition 
fluctuations, which will xertainly influence the performance and reproducibility in the fabrication of such devices. 
Another suggested alloyjj Sii_a;Ca;, presents a much larger lattice mismatch within the alloy layer, leading to a much 
stronger randomness in local strain, which implies a higher probability of defects and smaller critical thickness if one 
attempts to grow strained Si on top of such a substrate. 

We have demonstrated that it is possible to achieve hydrogenic behavior for substitutional donors positioufid exactly 
in the same (001) lattice plane parallel to the interface with the alloy substrate. For the recently proposedtEl "bottom- 
up" approach for Si QC fabrication, in which MBE growth is followed by positioning individual P impurity on the 
surface with the help of an STM tip, our result implies that small displacements of the P atoms due to surface diffusion 
or hops among the Si atoms forming dimers at the Si(001)2 x 1 surface should be acceptable for the rehability of 
QC operations mediated by exchange coupling. The flnal stage of the "bottom-up" approach involves Si overgrowth, 
encapsulating the deposited P donors. This is a crucial step for either relaxed or strained Si, since any atomic-scale 
^-component added to the relative position R leads to oscillatory behavior and possibly strong reduction of J. In 
MBE growth, surface binding and elastic energy differences among the various species lead to segregation effects, 
driving P atoms from previously grown layers into newly deposited ones. Segregation may be reduced by lowering 
the growth temperature and by increasing the growth rate,Eil but it may not be eliminated entirely. In a different, 
"top-down", approach, the P array in Si is formed by ion-implantation techniques, which implies a low degree of 
control over Rz on a lattice-parameter length scale. Residual variations in the relative positions of the donors must 
be corrected by individually calibrating the A and J gates. In other words, one must rely on probing and "correcting" 
each donor bound state through the surface gates. In general, our work establishes that Si QC design could tolerate 
some imprecision in the positioning of the dopant atoms in a 2D layer parallel to the interface, but any lack of precise 
control in the z direction could be fatal. 

We point out an additional advantage of strained over relaxed Si according to our numerical results. Following all 
overgrowth stages, one A gate must be placed on top of each buried impurity and one J gate in-between neighboring 
impurity pairs (see Fig. |l](a)). We have shown that the exchange coupling is enhanced for donor pairs parallel to the 
interface in negatively strained Si conflgurations. This allows for larger inter-donor distances as compared to relaxed 
Si, which is particularly convenient in terms of fabricating and accurately positioning the external gates. 

Finally, we point out the obvious fact that, although we have emphasized the relevance of our theory to Si QC 
architectures throughout this paper, our work is of general validity completely transcending QC considerations and 
context. In particular, we have developed a theory for calculating the donor exchange energy in strained Si systems 
incorporating the quantum interference effect arising from the multiple valley structure of Si conduction band. Our 
results should apply to any problem involving Si donor exchange coupling considerations. 
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APPENDIX A: HEITLER-LONDON APPROXIMATION FOR DONORS IN SI 

In this paper the Heitler-London approximation is employed to calculate the exchange splitting J between the 
ground singlet and triplet states for two ^^P donors embedded in strained Si. Here we derive the expression for J in 
terms of the donor relative position and the host band structure parameters. Assuming the ground donor electron 
wave functions to be ipir), the singlet and triplet wave functions for donors located at the origin and R are 
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Mri,r2) = ^^^^^ [V^CrOV-Cr^ - R) + ^(ri - RMr,)] x i= (| ti> - | it)) 
Mri,r2) - [^(ri)V(r2 - R) - ^{ri - R)V'(r2)] x i= (| ti) + | it)) 



where S" is the overlap integral 

The Hamiltonian that governs this two-electron system is 



e 
e 



|r. - R| 



ri-r2 R 



(Al) 



(A2) 



The energy splitting between the ground singlet and triplet states within the Heitler-London approximation id 

{HqS'^ - Hi) . 

The terms appearing both here and in dTh are 



Hi^ J d3rid3r2V*(ri)V*(r2 - R)0^(ri - R)V'(r2) 



with 



o = — 

e 



In-RI Iral |ri-r2| R 



(A3) 
(A4) 

(A5) 



Replacing the ground state wave functions of the form for in the expressions above leads to sums over integrals 
which involve different envelope functions and phases coming from the Bloch wave functions. For example, for the 
overlap (Al) one gets 

^M^KK' /d^ri^^^(ri)F.(ri-R)e-(''''+K)--^e^('^"+K')-("--^) , 



J dhiF^{ri)F,{ri - R)e^(k.-k,+K'-K).ri 



(A6) 



where 



A, 



i/KK' = a ai, Cj^ 



Following the approximations suggested in Ref. we neglect all integrals which contain rapidly oscillating phases in 
the integrated variables ri or r2. The remaining integrals involve different anisotropic envelope functions, and may be 
cast into forms similar to those appearing in the H2 molecule HL problem through a coordinate rescaling along the 
three Cartesian directions by the respective effective Bohr radii. For the specific example above, we get: 



IJ.K. 



(A7) 



with 



^ = e-^.(l + 7^^+7^^/3) 



(A8) 
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Equation ( A8) is exactly the hydrogenic expression for the overlap, but e valu ated a t th e rescaled internuclear distance 
which, for /i = z, is TZz = {Rxla,Ry/a,Rz/b). Integrals appearing in (A3) and (A4) contain a distance-dependent 
denominator which, at large internuclear separations, may be replaced by the value appropriate for the line joining 
the two donors .tfl For these terms it is convenient to define 
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1/2 



(A9) 



where 9^ is the polar angle between the /z-direction and R. 

Taking all these considerations into account, we arrive at the final expression for the exchange coupling J(R): 



^(R) = E 



E 

K,K' 



The expression for the kernel is 
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|2gi(K-K')-R 



^Jf,,{Il) cos(k^-k,)-R. 
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A 

C2=El"^l'l«7lM, 
A7 



(All) 

(A12) 
(A13) 



ca = /(0Al[-l/7^A 
/(^A.)[- 



-2K, 



(1 + 1/^a)], 



rAA 



/(^A)[l/7^A 



-2-Rx 



(l/7^A + 11/8 + (3/4)7^A + (l/6)7^2)], v^, 



-e ■'^^'(1 + 7^^)], s^, is given in ([A^), and j^^ = /(^A')i^f C^p)' '^it^^ J^f(^) given in (B5) of Ref. As in 
Ref. we make the following additional assumption regarding the off-diagonal elements of c' and y. = y'j^^jyy 

and c'i 



-A7 — V^aaSt 



The expression (AlO) for exchange J(R) is formally equivalent to the J(R) obtained by Andres et aZ0 under the 



additional approximation 4>fi(v) 

j^,(R) = y d3rid3i.2j7^*(i.i)^^;(r2-R) 



<^K,Oi but with J'^i, replaced by 
2 



e|ri - r2 



■F^(ri -R)^^,(r2). 



(A14) 



This result emerges from an approximation for the exchange coupling commonly adopted in the liter ature,aili 
which only the exchange integral of the Coulomb electron-electron repulsion term is included. To our knowledge, no 
formal justification for this assignment is available for electrons bound to donor impurities via Coulomb interaction. 
However, results presented in Fig. |^ show that this is actually an acceptable approximation if the inter-donor separation 
is between 100 and 200 A in Si, where relative changes |AJ(R)|/J(R) are at most 20%, keeping the same qualitative 
features and trends reporte d her e for all strain and R values investigate d. In fact, it is possible that the approximate 
expression obtained using ( A14 ), which is positive-definite, instead of ( All ), becomes more reliable than HL in the 
R ^ 00 limit, given the artifact of the HL expression for the exchange splitting in H2 molecule that predicts an 
unphysical negative J for R beyond 50 Bohr radii.oa 
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(a) Donor nuclear spin (b) Donor electron spin 

Si quantum computer SiGe quantum computer 




FIG. 1. (color) (a) Basic ingredients of the Silicon QC.H An array of substitutional ^^P atoms in Si is schematically 
represented in the figure by only two atoms. The active elements in the QC are the the electronic and nuclear spins in the ''^P 
donors. The nuclear spins (7 = 1/2) are the qubits. The electrons mediate an effective exchange interaction between nuclear 
spins and also participate in qubit initialization and readout operations. All operations are controlled by the A-gates, placed 
above each donor, and by the J-gates, placed between neighboring donor pairs. A barrier layer separates the host Si crystal 
from the metallic gates. The A-gate is biased to bring the corresponding nuclear spin in and out of resonance with an external 
rf magnetic field, allowing one-qubit logic operations to be performed. Donor electrons mediate the nuclear spin interactions 
through their exchange coupling. Electron exchange and consequently nuclear spin interactions are turned on orpoff by the 
J-gate. A uniform magnetic field ~ 1 — 2 Tesla is applied along the z direction, (b) Basic ingredients of the Si/Ge QC.u An array 
of substitutional ^^P atoms in an alloy layer of a Si/Ge composition-modulated heterostructure is schematically represented in 
the figure by only two atoms. The qubits are the donor electron spins. Given the difference in the electronic g-factor in Ge 
(g = 1.5) and in Si {g = 2.0), electron-spin resonance may be controlled by the top gate, which drives the corresponding donor 
electron into regions of different alloy compositions. Two-qubit exchange interactions are turned on through the top gates by 
drawing neighboring donor electrons into a layer of small effective mass and larger wavefunction overlap. 




X (valley strain) 



FIG. 2. (a) Energy splitting between the ground and first excited state of a single electron bound to a single ^^P donor in 
strained Si host as a function of the valley strain x- For Si commensurately grown on Sio.8Geo.2, with x ~ —20, this energy 
splitting is about 3.3 meV, down from ~ 12 meV in relaxed bulk Si. For x > 0, the first excited state is a doublet. At large 
X, this energy splitting is about 7.5 meV. (b)Hyperfine coupling of the ground single donor electron state, normalized to the 
unstrained host value, as a function of the valley strain x- At x = —20, the ground state hyperfine coupling is about 38% of 
the ground state in relaxed bulk Si. 
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Inter-donor distance (Angstrom) 

FIG. 3. Calculated donor electron exchange in Si as a function of inter-donor distance and valley strain x, with the two 
donors located along the [100] axis (with the same y and z coordinates). The symbols represent exchange values when both 
donors are on substitutional lattice sites, while the curves give results for continuously varied inter-donor distance. The filled 
symbol in particular refers to the configuration in which one donor is at the origin while the second donor is 20 lattice spacings 
away along x axis. This configuration will be used as a starting reference point for the next two figures. Notice that for systems 
with different values of x, the symbols representing lattice sites are slightly displaced horizontally relative to each other. The 
reason is self-evident, as the strain is produced by growing Si commensurately on substrates ranging from Si/Ge (x < 0) to 
Si/C (x > 0) alloys, with varying alloy compositions that lead to different lattice spacings. 
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FIG. 4. Calculated exchange coupling as a function of displacements of one of the donors relative to the reference config- 
uration in Si lattice uniaxially strained along z axis with x = —20. The origin represents the situation when the two donors 
are in the reference configuration defined in Fig. ^ The two curves in the figure refer to the cases when one of the donors 
is displaced from its reference site along the y and z axis, respectively. The symbols represent the substitutional lattice sites 
along the two axis. 
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FIG. 5. Calculated exchange coupling as a function of displacements of one of the donors relative to the reference config- 
uration in Si lattice uniaxially strained along z axis with x = —20. The origin represents the situation when the two donors 
are in the reference configuration defined in Fig. ^ The three curves in the figure refer to the cases when one of the donors is 
displaced from its reference site along the [110], [Oil], and [101] axis, respectively. The symbols represent the substitutional 
lattice sites along the three axis. The curve for the [Oil] axis is symmetric around the reference site because the corresponding 
displacement along the [Oil] axis causes symmetric variations in the inter-donor distance |R|. 




100 120 140 160 180 

Inter-donor distance (Angstrom) 



FIG. 6. Calculated exchange coupling for in-plane displacements of the donors in the x-y plane for relaxed Si (x = 0, 
diamond symbols) and strained Si (uniaxial strain along the z direction, x = —20, square symbols). The inset, in which the 
filled dots represent the donors, describes schematically the positions considered for the displaced donor, which consist of all 
possible lattice sites between two concentric circles of radii 90 A and 180 A with the other donor positioned at the center of 
the circles. The data points correspond to the exchange calculated at all relative positions considered. The solid line is J(R) 
for R along the [100] direction for x = —20. 
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FIG. 7. Comparison of the results from calculations with the Heitler-London approximation (solid curve) and its approximate 
form (dashed curve) for the donor exchange in bulk relaxed Si. Without loss of generality for the comparison, the two donors 
here are placed along [100] axis. The inset gives the difference between these two approaches in units of the HL value for 
exchange as a function of the inter-donor distance. At about 180 A donor separation there is a sign change in the exchange 
difference AJ, while in the range of distances of practical interest, from 100 to 200 A, the relative difference remains small. 
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